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We investigate the connection between a formal property of the critical behavior of several systems 
in the presence of quenched disorder, known as "dimensional reduction" , and the presence in the 
same systems at zero temperature of collective events known as "avalanches" . Avalanches generically 
produce nonanalyticities in the functional dependence of the cumulants of the renormalized disorder. 
We show that this leads to a breakdown of the dimensional reduction predictions if and only if the 
fractal dimension characterizing the scaling properties of the avalanches is exactly equal to the 
difference between the dimension of space and the scaling dimension of the primary field, e.g. the 
magnetization in a random field model. This is proven by combining scaling theory and functional 
renormalization group. We therefore clarify the puzzle of why dimensional reduction remains valid in 
random field systems above a nontrivial dimension (but fails below), always applies to the statistics 
of branched polymer and is always wrong in elastic models of interfaces in a random environment. 
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In the theory of disordered systems, "dimensional re- 
duction" is the property shared by some models that the 
long-distance physics in the presence of quenched disor- 
der in some spatial dimension d is the same as that of 
the pure model with no disorder in a reduced spatial di- 
mension, d — 2. In the known examples where it has been 
found through perturbation theory, i. e. the random field 
Ising model (RFIM)i~— , elastic manifolds in a random 
environment, the random field and random anisotropy 
O(N) models^— £ and the statistics of dilute branched 
polymers 5 -^, it entails two conditions: (1) that the long- 
distance physics is controlled by a zero-temperature fixed 
point, so that it can be equally described from the solu- 
tion^) of a stochastic field equation at zero temperature, 
and (2) that an underlying supersymmetry (involving ro- 
tational invariance in superspace) emerges in the field- 
theoretical treatment of the stochastic equation 5 -*^. Di- 
mensional reduction however is known to be wrong in 
some of the above models, the random field and random 
anisotropy models in low enough dimension (a rigorous 
proof exists for the RFIM in d = 2 and 32*i£) and the 
elastic manifold in a random environment^. On the other 
hand, it is proven to be right for the branched polymer 
case in all dimensions below the upper critical on o 11 ' 12 . 

We have recently shown that the breakdown of di- 
mensional reduction and the spontaneous breaking of the 
underlying supersymmetry take place below a nontrivial 
critical dimension-- in the random-field Ising, and more 
generally O(N), model: this dimension is close to 5 for 
the Ising (N = 1) version and decreases continuously as 
N increases until it reaches 4 when N approaches 18 (the 
upper critical dimension is equal to 6 for random field sys- 
tems . Describing this phenomenon requires a renor- 
malization group (RG) approach that is functional, as 
the origin of the dimensional reduction breakdown is the 
appearance of a nonanalytic dependence of the renormal- 
ized cumulants of the random field (a linear "cusp") in 
the dimensionless fields, and nonperturbative, as it takes 
place away from regimes where some form of perturbation 



analysis is possible (except for the O(N) model when d is 
close to the lower critical dimension of l 13 ' 15 ) . A similar 
conclusion was previously reached for the case of an elas- 
tic manifold in a random environment, but in this model 
the dimensional reduction predictions fail for all dimen- 
sions at and below the upper critical dimension (here 
equal to 4) and can be already assessed through a func- 
tional but perturbative RG^— . What is the physical 
mechanism behind these seemingly formal results? 

The existence of a cusp in the cumulants of the renor- 
malized disorder can be assigned to the presence of 
collective events known as "avalanches" . In any typ- 
ical sample of a disordered model, the ground state, 
which is the relevant configuration that describes the 
equilibrium properties of the system at zero tempera- 
ture, abruptly changes for specific values of the external 
source; the location of these abrupt changes are sample- 
dependent and the configurational change between two 
ground states is precisely an avalanche^—. (The latter 
is sometimes called a "static" avalanche^ 9 -.) The same 
phenomenon is observed, still at zero temperature, when 
the system is driven by the external source without be- 
ing allowed to equilibrate. The corresponding "dynamic" 
avalanches then take place out of equilibrium, between 
two metastable states of the syste m 25 ' 29 " —. 

However, the fact that abrupt changes corresponding 
to discontinuous variations of the magnetization (to use 
the language of magnetic systems) are found at zero tem- 
perature should come as no surprise. In disordered sys- 
tems, this can take place even in noninteracting zero- 
dimensional models. Consider for instance a d = (sin- 
gle point) c/) 4 theory with parameters such that the po- 
tential has two minima and couple the field to a ran- 
dom source h and a controllable source J. Then, ac- 
cording to the value of h + J, the ground state of the 
system will switch from the vicinity of one minimum to 
that of the other one with a jump in the magnetization 
(the field 4>) when h + J = 0. This jump corresponds 
to an avalanche, albeit a zero-dimensional one. As will 
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be illustrated in more detail below, these avalanches do 
generate a cusp in the cumulants of the effective random 
field (see Fig. 1). However, avalanches, and the result- 
ing nonanalyticities in renormalized disorder cumulants, 
affect the long-distance physics of a d-dimensional disor- 
dered model and play a role in the failure of dimensional 
reduction only if they are of collective origin and occur 
on sufficiently large scales. 

From the above discussion we conclude that in dis- 
ordered systems, (i) the physics of dimensional reduc- 
tion and its breakdown is associated with an effective 
zero-temperature theory, (ii) avalanches corresponding 
to discontinuous jumps of the magnetization and more 
generally of the field configuration are the rule rather 
than the exception at zero temperature and (iii) the 
dimensional-reduction breakdown results from the pres- 
ence of avalanches that must be collective enough to influ- 
ence the long-distance properties of the model. The cen- 
tral question is then: Under which conditions does this 
happen? We show that dimensional reduction remains 
valid when the exponent df that characterizes the scal- 
ing behavior of the largest typical avalanches at criticality 
is strictly less than the difference between the spatial di- 
mension d and the scaling dimension of the field eL near 
the relevant zero-temperature fixed point: df < d — d$. 
This condition is satisfied for the RFIM at and close to 
the upper critical dimension d uc = 6 and for branched 
polymers at and around the upper critical dimension 
d uc = 8 and below d = 4 (at least). On the other hand, 
dimensional reduction breakdown takes place when df is 
equal to d — d^r°. This is found for elastic manifolds 
in a random environment at and near the upper criti- 
cal dimension d u = 4 as well as for the RFIM below a 
critical dimension close to 5. This clarifies the intriguing 
result of why dimensional reduction fails below a nontriv- 
ial dimension for the RFIM, is always broken for random 
elastic manifolds, but applies to the branched-polymer 
problem. 

Avalanches and their consequences on the cu- 
mulants of the renormalized disorder. Take a dis- 
ordered system at zero temperature in which avalanches, 
representing discontinuous changes of the relevant con- 
figuration of the system under a variation of an exter- 
nal source, are present. We use the language of mag- 
netic systems and characterize configurations by the lo- 
cal magnetization. All considerations, however, equally 
well apply to configurations described by a continuous 
field, in or out of equilibrium, and to nonmagnetic sys- 
tems, e.g. to the displacement of an interface in a random 
mediu m 20 i 21 i 29 i 35 or to the appropriate density field de- 
scribing the universal properties of dilute branched poly- 
mers. We focus on situations in which the field under 
study (the "magnetization" ) is the local order parameter 
and is linearly coupled to the external source, but one 
should keep in mind that avalanches could also be trig- 
gered by changing for instance the strength of the dis- 
orde r 36 ' 37 at zero temperature and that avalanches may 
involve a field that is not the local order parameter (e.g., 



the magnetization in a spin glas o 38-40 ). 

Consider for simplicity an external source J that is 
uniform in space. The avalanches can then be charac- 
terized by their size S (the overall change in the total 
magnetization^) whose distribution is described by a 
density p(S,J), such that p(S, J)dS dJ is the (disorder 
averaged) number of avalanches of size between S and 
S + dS when the source is between J and J + dJ. The 
magnetization m( J; h) is the spatial average of the local 
order parameter field for a given sample characterized by 
the disorder realization h. Its change between two values 
of the external source J\ and J 2 is the sum of two con- 
tributions: a first one comes from the smooth changes 
in the ground state (in equilibrium) or in the metastable 
state (out of equilibrium) and another one comes from 
the avalanches that take place between J\ and J 2 (with 
Ji> J\)- As a consequence, the moments of the differ- 
ence [m(Ji;h) — m(J2;h)], which is a random variable, 
are given by 

[m(J 2 ;h) -m(Ji;h)]P = 

1 f j2 f 00 (!) 
(j^T J j dJ'J^ dSSVp{S,J')+reg, 

where the first term is due to avalanches taking place 
at the same value of the external source and the second 
term, denoted by reg, includes the contributions involv- 
ing the smooth variation of the magnetization or distinct 
avalanches. S m i n is a microscopic lower cutoff on the 
size of the avalanches, L d is the sample volume, and the 
overline denotes the average over the quenched disorder. 

For even moments, due to the symmetry in the ex- 
change of Ji and J 2 , the first term in Eq. ([I]) gives rise 
to a linear cusp when J2 — > Jv 

[m(J 2 ;h)-m(J i; h)] 2 P 

= I J 2 - Ji\L- 2 P d / dS S 2p p(S, J) + 0([J 2 - J1] 2 ). 

Js min 

(2) 

[The other term in Eq. (TTJ) is regular or at least less sin- 
gular than the first one and is therefore a 0([J 2 — >/i] 2 )-] 
It is easily realized that the pth moment is obtained by 
considering disorder averages over p copies of the same 
sample, with each copy coupled to a distinct external 
source J a , a = 1, • ■ • , p. The first moment is the change 
in the average magnetization of the system and the pth 
cumulant (which is the connected piece of the pth mo- 
ment) is then related to Green's functions at zero mo- 
mentum of a p-copy system in which each copy is char- 
acterized by a distinct source. For example, the relevant 
2-point Green's function, given by G(q — Q\J\,J 2 ) — 
L d [m(Ji; h)m(J2; h) — m(Ji;h) m(J 2 ;h)], is an exten- 
sion to general sources J\ ^ J2 of what is usually called 
the "disconnected" 2-point correlation function in the 
theory of disordered systems. One-particle irreducible 
(1PI) correlation functions (or proper vertices)^ associ- 
ated with the above Green's functions can be introduced 
along the same finest. From Eq. @, one immediately 
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derives that, for instance, the 2-point Green's function 
G(q = 0; J\ = J — 8 J, J 2 = J + 8 J) has a nonanalytic 
dependence as 8 J — > 0, with the amplitude of the cusp 
related to the second moment of the avalanches: 



G(0; J-8J,J + 8 J) - G(0; J, J) = 
1 f°° 

-\8J\j^J dS S 2 p(S,J) + 0{SJ 2 ). 



(3) 



This can be transposed to the associated 1PI vertices and 
can generalized to higher orders as well 6 -^. 
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FIG. 1: Illustration of avalanches and their conse- 
quence on the functional dependence of the Green's func- 
tions in the schematic case of the d = RFIM stud- 
ied in equilibrium at T = 0. (a) Free energy F(<fi) — 
(J + h)<f> versus <f> for different values of J, with T(4>) — 
— ( |t |/2)0 2 + (g/4!)</> 4 ; (b) Ground state configuration associ- 
ated with (a); Two-point Green's function G(0; — 8 J, 8 J) = 
[4>as(-SJ + h)</>Gs(SJ + h) - (f>Gs(-SJ + h) cj> as (8J + h)}, 
where the average is over a Gaussian distributed random field 
h: notice the linear cusp around 8 J = 0. 

The above arguments show that avalanches induce a 
linear cusp in the functional dependence of the correla- 
tion functions describing the cumulants of the effective 
or renormalized disorder at T — 0. This is always true 
even when the avalanches take place on a restricted scale, 
away from any critical conditions (see Fig. 1 for an illus- 
tration). However, we are interested in the long-distance 
behavior of disordered systems. We therefore study the 
situation in which avalanches occur on all scales, as found 
for instance in the RFIM at the equilibrium or out-of- 
equilibrium critical point, in the rough phase of an elastic 



manifold pinned by a random medium or at its depinning 
transition, etc. 

Consider a d-dimensional system of linear size L at 
zero temperature. At large scale, when the correlation 
length and the extent of the largest typical avalanches 
have reached the system size, one expects that the (nor- 
malized) probability density of having an avalanche of 
size S (see footnote 5) can be written in the following 
scaling form^i^ 6 -^^: 

D L (S,J)=S- r V(-f,\J-J e \S*) (4) 

where Sl ~ L df is the size the largest typical "crit- 
ical" avalanches in the finite system 64 ; Sl acts as a 
cutoff for the scaling function T> that decays exponen- 
tially for S/Sl > I- Critical conditions correspond to 
J = J c (for the RFIM at equilibrium one has J c = 
due to the Z 2 symmetry and for the random manifolds 
there is no condition on J as the whole phase is crit- 
ical). The avalanche size distribution is normalized so 
that J s dS Dl(S, J) — 1 and the moments of the nor- 
malized avalanche size distribution are then defined as 
< SP > L = J™ min dS S p D L (S,J). When 1 < r < 2, 
which is usually found (for instance, the mean-field value 
for r is equal to 3/2 in all models), the normalization 
factor is dominated by the small avalanches whereas all 
moments with p > 1 are dominated by the largest typical 
avalanches and behave as (Sl) p+1 ~ t ~ L^ p+1 ^ T ^ df when 
L — > oo. 

We keep using the language of magnetic systems and 
let ?til( J; h) denote the magnetization of a given sample 
of linear size L. (Here and below, we explicitly indicate 
the dependence on the system size L\ this makes the 
expressions somewhat clumsy but will be helpful later 
on to make the connection with the results of the func- 
tional RG.) We are primarily interested in the moments 
of the random variable nii( J 2 ] h) — mt{J\ \ h). The den- 
sity of avalanches [see Eq. ((T|)] is related to the nor- 
malized probability density by an overall L and J de- 
pendent factor: pl(S,J) = Po,l(J)Dl(S, J). As dis- 
cussed before in connection to Eq. ([1]), the first moment 
of mL(J2jh) — mi(Ji;h) contains a contribution from 
the smooth change of the magnetization and one from 
the avalanches. The so-called "connected" susceptibility 
Xc,l{ J), which is the standard magnetic susceptibility di- 
vided by the temperature in order to have a proper zero- 
temperature limit and which is obtained by deriving uil 
with respect to J, can then be expressed as 

/ t\ smooth/ j\ , 1 

XcA J ) = Xc, L {J) + Yd 



dSS PL (S,J). (5) 



Under critical conditions, Xc.l goes as L 2_I? . By fur- 
ther making the natural assumption that the contribu- 
tion from the avalanches is of the order or larger than 
the smooth one and by using Eqs. (|4I5[) as well as the 
fact that the first moment of the avalanches < S >l is 
dominated by large avalanches, one then obtains that 

Po, L (Jc)L- d+i2 - T)df ~ L 2 -\ (6) 
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As a result, Eq. ([2]) can be rewritten as 

[m L (J 2 ;h)-m L (Ji;h)]P 

~ | J 2 - Ji| tf-i-to-W-W + 0((J 2 - Ji) 2 ) 

and the linear cusp in the 2-point Green's function at 
zero momentum Gl (q = 0; J — 6 J, J + 8 J) when 8 J — > 
[see eq. ([3j] is found as 

Gl(9 = 0; J c - 8 J, J c + 8 J) -G L {q = 0; J c , J c ) 
~ |JJ| i d /+ 2 -" ( j 

up to a 0(8J 2 ). The amplitude of the cusp therefore di- 
verges as the size of the system diverges. (One should 
however keep in mind that the whole function Gl{q = 0) 
itself diverges as L 4 ~ r ' at criticality.) The above re- 
sult generalizes to higher-order Green's functions through 
their relation to the cumulants of the magnetization. 

From Eq. ([5]), it is easily derived that the associated 
1PI correlation function A^(q = 0; mi, OT2), which is the 
second cumulant of the renormalized disorder, also has 
a cusp in \rri2 — mi\ as m 2 — > mi. After introducing 
mi = m c — 8m and m 2 — m c +8m, where 8m — > and m c 
corresponds to the value of the average magnetization at 
criticality, and using the relation between Green's func- 
tions and 1PI functions (see ref4i and footnote 6) as well 
as 8m = mh{Jc + 8 J; h) - tol(J c - 8 J; h) ~ SJxc,l(J c ) 
when 8 J — > 0, we obtain, up to a 0(8m 2 ), 

A L (q — 0; m c - 8m, m c + 8m) - A L (q = 0; m c , m c ) ~ 
XcAJcr 2 [G L {q = 0; J c - SJ, J c + 8 J) — Gh{q = 0; J c , J c )\ 
~ \Sm\ L d f- 2 ^ 

(9) 

where we have also used that Xc,l(J c ) ~ L 2 v - Again, 
the above expression can be generalized to higher-order 
cumulants of the renormalized disorder, which all display 
a cusp in their functional dependence with an amplitude 
that is system-size dependent at criticality. 

Functional RG and dimensional reduction 
breakdown. As already stressed, the functional RG is a 
powerful and necessary framework to describe the criti- 
cal behavior of the disordered systems of interest. Within 
such an approach, which is a version of Wilson's continu- 
ous RG^2r— , the fluctuations are progressively taken into 
account by introducing an infrared cutoff that enforces 
the decoupling of the low- and high-momentum modes 
at a running scale k. For k = 0, all fluctuations are in- 
cluded and the exact theory is recovered. One ends up 
with flow equations for the moments of the renormalized 
disorder that describe the evolution of these moments as 
one decreases the infrared scale k. For instance, an equa- 
tion is obtained for the second cumulant of the renormal- 
ized random field or random force A k (q = 0; (j>\, 2 )^~— , 
which is the quantity already considered in the previous 
sections, with the field 4> and the infrared scale k playing 
here the same role as the local magnetization m and the 
inverse system size 1/L, respectively. (This RG equation 



for Afc is in general part of a hierarchy of coupled flow 
equations.) 

In order to reach the fixed point that controls the long- 
distance behavior under study, one must introduce scal- 
ing dimensions and convert the quantities appearing in 
the RG flow equations from "dimensionful" to "dimen- 
sionless" . For the cases of interest where on the one hand 
avalanches are present and on the other hand dimensional 
reduction is found in standard perturbation theory, we 
have stressed that the fixed point is at zero temperature. 
Temperature is then a dangerously irrelevant variable, 
and an associated exponent 6 > is introduced through 
an appropriate definition of a renormalized temperature 
T«i2ri2i£^: T k ~ k e '. Near the zero-temperature fixed 
point, the dimension of the field cj> is modified from its 
standard value of (d — 2 + 77) /2, with r\ the anomalous di- 
mension, by a term involving the temperature exponent: 

dt = ^(d-2 + v -e) = ±(d-4 + fj), (fO) 

where we have also introduced the additional anomalous 
dimension fj through the relation 9 = 2 + 77 — 77. Similarly, 
the second cumulant of the renormalized random field or 
force has the scaling dimension of a 2-point IPI vertex, 
2 — 7/, modified by the temperature exponent, i.e. 2 — 77 — 
6 = —2ij + fj; it can be put in dimensionless form as 

A k (q = 0; cf>i, 0a) - fc- (2 "-^ fc (0; <pi, ^ 2 ), (f f ) 

where tp is the dimensionless field (see also SI appendix) . 

Dimensional reduction corresponds to 9 = 2, which 
implies fj — m and to all other exponents equal to their 
value in the system without disorder in dimension d — 
2. The main outcome of the functional RG studies is 
that breakdown of dimensional reduction is related to 
the presence of a cusp in the functional dependence of 
the dimensionless second cumulant of the renormalized 
random field or force, <5fc(0; ipi, (P2), in the vicinity of the 
zero-temperature fixed pointer—. More concretely, after 
introducing ip = (ipi + <^2)/2 and Sep = (tp2 — ¥>i)/2, 
the "cuspy" behavior that changes the critical exponents 
from their dimensional reduction prediction is of the form 

5,(0; ip - Sip, ip + Sip) = 6*fl(tp) + 8* >a {ip)\8ip\ + 0(8tp 2 ) 

(12) 

when Sip — > 0, with (5». a < 0; the star indicates the fixed- 
point value at k = 0. As a result of a nonzero <5*. a , the 
exponent 9 takes a nontrivial (^-dependent value < 2 and 
77 and 77 differ from the dimensional reduction values, with 
77 ^ 77. 

The connection between the quantities computed 
through the functional RG and those discussed in the 
previous sections can be made by associating the infrared 
cutoff k with the inverse of the linear extent of the sys- 
tem, i.e. k ~ 1/L. Eq. ([9]) can then be expressed in a 
dimensionless form by dividing the cumulant and the 
magnetization 5m by their scaling dimensions L 2ri ^ ri and 
L~~ ( d_4 +'5)/ 2 , respectively. We immediately obtain that 
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the amplitude of the linear cusp in dimensionless form 
scales as with 

L d s -2(2-7 1 )-(2r,-f))-^(d-i+n) _ L d f -±(d+4-fj) ^ / 13 \ 

which can also be rewritten as L df ~( d ~ d,f '\ By compar- 
ison with Eq. (|12[) . one therefore finds that the cusp 
persists in the dimensionless quantities when L — > oo, 
i.e. at the fixed point, if and only if df = d — d<p. If 
df < d — d(j>, the cusp is only subdominant and does not 
affect the leading critical behavior and the associated ex- 
ponents. (Note that the condition df > d — d$ is not 
compatible with the result of the functional RG studies, 
in which proper renormalized theories have always been 
found with no stronger nonanalyticities than the linear 
cusp.) 

We conclude from the above derivation that di- 
mensional reduction breaks down due to avalanches if 
the fractal dimension of the largest typical "critical" 
avalanches satisfies the condition 

d f = d-d$. (14) 

On the other hand, dimensional reduction remains valid 
if 

df < d — d^, (15) 

despite the presence of the avalanches and of a cusp in the 
dimensionful cumulants of the effective disorder at zero 
temperature. In the latter case, the difference (d — d^, — 
df) can be reinterpreted and computed in the functional 
RG. Indeed when perturbing the "cuspless" fixed-point 
valued of the dimensionless cumulant with a function 
that itself displays a linear cusp, the amplitude of the 
cuspy perturbation should go to zero as k —> in such a 
way that 

Sk{0; <p-dip, ip + 5(p) ~ 5*(0;(p — 5<p, ip + 5ip) + k x f\(ip, Sip) 

(16) 

with A = (d — d^) — df > and f\(ip,5tp) ~ \6ip\f\(ip) 
when dip — > 0. Information on the fractal dimension of 
the largest typical avalanches at criticality can then be 
obtained from an investigation of the irrelevant direc- 
tions associated with nonanalytic eigenfunctions around 
the fixed point. This is what we have done by solving 
the nonperturbative functional RG equations derived in 
RefJ£ for a function 5k of the form given in Eq. (fT6|) (see 
SI appendix). 

Results and discussion. We are now in a position 
to discuss the consequences of Eqs. (fl"4")l and (fTS"]) for sev- 
eral disordered systems in which dimensional reduction 
is predicted by standard perturbation theory. 

Consider first the mean- field limit. Avalanches are 
present at T = and the distribution of the avalanche 
sizes at criticality can be described by the scaling ex- 
pression in Eq. The exponents r and dt can be 
easily derived for fully connected models. This was first 
done by Dahmen and Sethna - for the out-of-equilibrium 



behavior of the slowly driven RFIM at zero tempera- 
ture, but can be generalized to other models as well. 
(As usual, the mean-field exponents are expected to be 
"super- universal" .) The values of the avalanche expo- 
nents are r = 3/2 and df = 4. These values have 
also been recovered by Le Doussal and Wiese^ from a 
field-theoretical treatment of elastic manifolds in a ran- 
dom environment. What conclusion can then be drawn 
about the influence of the avalanches on the long-distance 
physics? At the upper critical dimension d uc at which the 
exponents take their mean-field values, the anomalous di- 
mension f\ — so that d^ = (d uc — 4)/2. The conditions 
in Eqs. (TT4|) and (fTS)) simply amount to comparing df = 4 
and duel '2 + 2. For random field and random anisotropy 
models (we include here models with iV-componcnt fields 
which we expect to behave in a similar manner as that of 
the single-component one), the upper critical dimension 
is d uc — 6 so that dimensional reduction should apply. 
The same is true for the statistics of dilute branched poly- 
mer for which the upper critical dimension is d uc = 8 49 . 
On the other hand, for interfaces in a disordered environ- 
ment, the upper critical dimension is d uc = 4: a failure 
of dimensional reduction is then expected, possibly only 
in logarithmic corrections at d = d uc but more severe as 
one lowers the dimension. 

As one decreases the dimension from the upper criti- 
cal one, there must be a nonzero range of dimensions for 
which the dimensional reduction predictions correctly de- 
scribe the critical behavior of the random field, random 
anisotropy and branched polymer models, but likely not 
that of the random manifold one. Actually, the latter 
model has been studied in great detail by the perturba- 
tive functional RG in d = 4 — c 21 ' 29 i 35 . It was found 
that the dimension d / characterizing the cutoff on large 
avalanches in the presence of a finite-size or infrared cut- 
off on the system is equal to d + £, where ( is the ex- 
ponent describing the roughness of the interface (for a 
single-component displacement field). As the dimension 
d(f, of the field is itself equal to — C (fj is formally equal to 
4 — d — 2£) , it follows that the equality df = d — d^ is al- 
ways verified and that dimension reduction never applies, 
as indeed found by direct computation of the critical ex- 
ponents within the functional RG or in computer sim- 
ulations. This conclusion is valid for the pinned phase, 
in equilibrium, and for the depinning threshold in the 
driven case. 

For the RFIM at equilibrium, we have shown through a 
nonperturbative functional RG that dimensional reduc- 
tion breaks down below a nontrivial critical dimension 
d cusp ~ 5.1 (see Refs i 14 ' 15 and footnote 2). According to 
the above conditions, the avalanche exponent df should 
then be equal to d — d^ = (d + 4 — fj)/2 below d cusp and 
to (d + 4 — fj)/2 — A, where A is the eigenvalue associ- 
ated with the irrelevant cuspy directions around the cus- 
pless fixed point (see preceding section and SI appendix), 
above d cusp . In Fig. 2, we plot the theoretical prediction 
for df based on the above relations and the computa- 
tion of dtj, and A from the solution of the flow equations 
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previously derived in our nonperturbative functional RG 
approach of the RFIM 15 (see also footnote 2). The pre- 
diction is confirmed at the upper (see above) and lower 
critical dimensions. For the latter, di c = 2, one indeed 
expects the avalanches to be compact even at criticality 
and their fractal dimension therefore to be equal to the 
spatial dimension, dt — d — 2 (see also Ref^i). As the 
dimension of the field is equal to zero, d^ = (and fj = 2), 
the equality in Eq. (fT4|) is satisfied. (Note that the re- 
sults at the lower and upper critical dimensions apply to 
both the equilibrium and out-of-equilibrium critical be- 
havior of the RFIM 66 .) Beside this, direct measurements 
or computations of the avalanche exponent df are un- 
fortunately scarce. We therefore suggest that, as done 
for models of an elastic interface in a disordered environ- 
men t 20 ' 21 ' 29 and for the out-of-equilibrium, metastable 
behavior of the driven RFI M 25 ' 31 ~ — , systematic studies 
of the avalanches and of the cumulants of the effective 
(renormalized) disorder would be worthwhile, e.g. in the 
ground state of the RFIM, and would allow a direct test 
of the predictions made on the basis of the functional 
RG. 

Finally, for the statistics of dilute branched polymers, 
so long as dimensional reduction applies, fj = 2rj and is 
negative. In consequence, d — d^ = d + A — fj > d + A. As 
the fractal dimension df should also be less than the di- 
mension d of spaced, one can see that df < d < (<i+4)/2 
when d < 4. From the condition on the scaling of the 
avalanches in Eq. (|15|) . we therefore obtain that dimen- 
sional reduction applies, at least, when d < 4 and in 
the vicinity of the upper critical dimension d uc = 8 (see 
above) ; the existence of an intermediate range of dimen- 
sions characterized by dimensional reduction breakdown 
is highly unlikely. This is of course in agreement with 
the known results according to which dimensional reduc- 
tion (to the Lee- Yang edge singularity, or equivalently to 
the universal repulsive gas singularity, in two fewer di- 
mensions) is always valid in the case of dilute branched 
polymer o 11 ' 12 . 

To summarize, we have related the breakdown of "di- 
mensional reduction" to the scaling characteristics of 
"avalanches". The former is the formal property ac- 
cording to which the critical behavior in the presence 
of quenched disorder is the same as that of the clean sys- 
tem in two dimensions less and is found within conven- 
tional perturbation theory in models whose long-distance 
physics is controlled by a zero-temperature fixed point, 
whereas the latter are large-scale physical events taking 
place in the relevant configuration of the system (ground 
state in equilibrium or metastable state out of equilib- 
rium) under the variation of an external source. This 
provides a solution to the puzzle of why dimensional re- 
duction breaks down in some models and not in others, 
in some range of dimensions and not for others. 

Note finally that at small but nonzero temperature, 
there are no avalanches and the variation of the relevant 
configuration of the system under a change of the exter- 
nal source is continuous, except possibly in mean-field 




cusp 



FIG. 2: Fractal dimension df of the largest typical avalanches 
versus space dimension d for the RFIM at the equilibrium crit- 
ical point, as predicted by Eq. (|14[) (in the region "noDR") 
and Eq. (|15[) (in the region "DR") and the nonperturbative 
functional RG. The numerical resolution of the RG flow equa- 
tions becomes extremely difficult in low dimension, typically 
for d < 2.9 and when approaching don — d cusp ~ 5.15 so 
that we have no result for these ranges of d. The filled cir- 
cles indicate the known values at the lower and upper criti- 
cal dimensions. The crosses are the numerical estimates for 
the out-of-equilibrium critical behavior of the driven RFIM in 
d = 3 25 ' 31 ' 33 and the square that for the equilibrium behav- 
ioi>2£. The dashed line is the upper bound (df < d). 



models. The nonanalyticities in the functional depen- 
dence of the cumulants of the renormalized disorder and 
of the associated Green's functions are then rounded in 
"thermal boundary layers" 14 ' 51 " —. In systems at equi- 
librium, these boundary layers are linked to the pres- 
ence of low-energy excitations that may also take place 
on large scales at and near criticality and are described 
as "droplets" 54 ' 55 . The relation between droplets and 
avalanches in disordered systems is by itself a very in- 
teresting topic, which however we have not considered 
here. In any case, this underscores that properties such as 
dimensional reduction and its breakdown or avalanches 
crucially depend on the system being at zero tempera- 
ture or having its critical behavior controlled by a zero- 
temperature fixed point. 

[Nonperturbative functional RG for the RFIM] 
We summarize here the main features of the nonper- 
turbative functional RG description of the equilibrium 
critical behavior of RFIM developed in Refs* 1 ^— as well 
as its extension to compute the stability of the zero- 
temperature fixed point against nonanalytic perturba- 
tions. 

The central quantity is the so-called "effective average 
action" Tj~i^ in which only fluctuations of modes with 
momentum larger than an infrared cutoff k are effectively 
taken into account. In the language of magnetic systems, 
r^, is the Gibbs free-energy functional of the local order 
parameter field obtained after a coarse-graining down to 
the (momentum) scale k. The effective average action 
obeys an exact RG equation under the variation of the 
infrared cutoff fc 4 ^. 
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In the presence of disorder, the generating or free- 
energy functionals are sample-dependent, i.e. random, 
and should therefore be characterized either by their 
probability distribution or by their cumulants. The 
latter description is more convenient as it focuses on 
quantities, cumulants and associated Green's functions, 
which are translationally invariant and can be generated 
through the introduction of copies (or "replicas" ) of the 
original system that are submitted to distinct external 
source o 13 ' 15 . The effective average action that generates 
the cumulants of the renormalized disorder, rfe[{<^> }], 
then depends on the local order parameter fields asso- 
ciated with the various copies a. It satisfies the following 
exact functional RG flow equatio n 13 ! 15 : 

d k T k [{fa}] = 

mx***>([*+*r)l- (17) 

(2) 

where is the matrix formed by the second functional 
derivatives of 1^ with respect to the fields 4> a (q) and 
Rf(q 2 ) = R k (q 2 )S ab + Rk(q 2 ), where R k and R k are 
infrared cutoff functions that enforce the decoupling of 
the low- and high-momentum modes at the scale k (see 
below). Note that the formalism can be upgraded to a su- 
perficld theory in order to describe the physics directly at 
zero temperature (i.e., at equilibrium, the ground-state 
properties)^ 5 -: this allows one to make the underlying su- 
persymmetry of the model explicit and to consider fea- 
tures associated with spontaneous or explicit breaking of 
the latteri 5 .. 

From Eq. (|17[) . one can derive a hierarchy of coupled 
RG flow equations for the cumulants of the renormalized 
disorder, r^i^i], T k 2[fa-, fa], etc, that are obtained from 
rfc[{0a}] through an expansion in increasing number of 
unrestricted sums over copies: 

Tfc [{fa}] = r ^ W -\Y, T ^[fa,fa] + ■■■ (18) 

a a.b 

At the microscopic scale, say k = A, the effective average 
action reduces to the "bare" action (or effective hamilto- 
nian) of the multy-copy system 

IW [{fa}] = I d d xY {\[dfa(x)] 2 + T -fa{x) 2 

+ j0a(*) 4 )-^ / d d xJ2<t>a(x)fa(x), 
' J a,b 

which generates the cumulants of the renormalized dis- 
order at the mean-field level (Ab is the bare variance of 
the random field that is taken with a Gaussian distribu- 
tion of zero mean). At the end of the flow, when k = 0, 
all fluctuations are incorporated and one recovers the ef- 
fective action r[{<^ Q }] which is the generating functional 
of the cumulants of the renormalized disorder and corre- 
sponds in the language of magnetic systems to the exact 
Gibbs free-energy functional of the multi-copy system. 



The detour via the superfield formalism provides a 
nonperturbative approximation scheme for the exact RG 
equation, Eq. (fT7|) . and a relation between the cutoff 
functions R k (q 2 ) and R k (q 2 ) which, both, do not explic- 
itly break the underlying supersymmetry at the origin 
of dimensional reduction. The minimal truncation of Y k 
that already contains the key features for a nonpertur- 
bative study of the long-distance physics of the RFIM is 
the following: 

r fc [{fa}} = J d d xJ2S^Z k (fa(x))[dfa(x)] 2 + U k (fa(x)) 
-\j d d x^V k (fa(x),fa(x)), 

(20) 

with three functions Z k , U k and V k to be determined. 
On the other hand, the cutoff functions must satisfy 
the relation R k (q 2 ) = -(A k J Z k )d q 2R k (q 2 ), with A fc 
the strength of the renormalized random field and Z k 
the field renormalization constant. Inserting the above 
ansatz for T k [{fa}] into Eq. (fT7|) leads to a set of cou- 
pled flow equations for the three functions Z k (fa), U k (fa 
and V k (fa,fa) (or alternatively its second derivative 
A k (fa, fa ) — d<p 1 d ( p 2 V k (fa,fa) which is the second cumu- 
lant of the renormalized random field at zero momentum 
discussed in the text). The RG is functional as its central 
objects are functions instead of coupling constants. 

One more step is needed to cast the nonperturbative 
functional RG flow equations in a form that is suitable 
for searching for the anticipated zero-temperature fixed 
points describing the critical behavior of the RFIM. One 
has to introduce appropriate scaling dimensions. This 
requires to define a renormalized temperature T k which 
should flow to zero as k — > 0. (This is the precise meaning 
of a "zero-temperature" fixed point.) Near such a fixed 
point, one has the following scaling dimensions: 

T k ~ k\ Z k ~ k-\ fa ~ k^ d - i+f >\ (21) 

with 9 and fj related through 9 = 2 + r\ — fj, as well as 

U k ~ k d - 9 , V k ~ k d - w , (22) 

so that the second cumulant of the renormalized random 
field A fc scales as k~^~^. 

Letting the dimensionless counterparts of U k , V k , A k , <p 
be denoted by lower-case letters, u k ,v k ,5 k , ip, the result- 
ing flow equations can be symbolically written as 

dtu' k ((f) = Pu'(<p), 

dtz k (<p) (23) 
d t S k ((pi,ip 2 ) = Ps(<Pi,<P2), 

where t = log(fc/A). The beta functions themselves de- 
pend on u' k , z k , S k and their derivatives [in addition, the 
running anomalous dimensions rj k and fj k are fixed by the 
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conditions Zk(0) = (5^(0,0) = 1]. Their expressions are 
given in RefJ£. 

Fixed points are studied by setting the left-hand sides 
of the equations in Eq. (|23lto zero. The zero-temperature 
fixed point controlling the critical behavior of the RFIM 
is once unstable and has been determined in a previous 
investigation^. We found that above a dimension close 
to 5.15, there exists a fixed point with no cusp singularity 
in the functional dependence of the associated 5*(<^i, ^2)- 
After introducing ip = (pi + p 2 )/2 and Sip = (ip2 — ipi)/2, 
the dimensionless second cumulant of the renormalized 
random field can indeed be expanded as 

S*(<p - Sp,ip + Sip) = <S*,o(v?) + ^*a{ { P)^ i P 2 + 0{\Sp\ 3 ) 

(24) 

when Sip — > 0. It can be shown that the critical behavior 
then satisfies dimensional reduction: fj — rj, so that 9 
1. and the critical exponents are exactly given by those 
of the pure Ising model in two dimensions less within the 
nonperturbative approximation that is the counterpart 
of Eq. dlJJl. 

To compute the eigenvalue A that characterizes the sta- 
bility of the above cusplcss fixed point with respect to a 
"cuspy" perturbation, we have considered the vicinity of 
the fixed point with u' k and Zk set at their fixed-point 
values and Sk(ip — Sip, ip + Sip) ~ 5*(ip — Sip, ip + Sip) + 
k x f\(tp,Sip) with fx(ip,Sip) ~ \Sip\fx{ip) when Sip -> 0. 
By linearizing the flow equation for Sk around <5* and ex- 
panding around Sip = it is easy to derive that f\(ip) 
satisfies the following eigenvalue equation: 

AAfoO = \iA - 4 + 3v)f\{<f) + \id - 4 + v)<Pf\{<p)+ 

f °° d ( 3 / 

v d d t dyy^ 1 l-f x (p)[4z'^p)p^(y,p)pf' 1) {y,p)+ 

4[z*(ip) + S '(y)} P ( ° A) (y, pf + [z'Kip) - S*, 2 (v)]p*(y, p) 2 ) + 
3fx(<p)p*(y, <p)(2[z*(tp) + s'(y)]pl 0,1) (?/, ip) + zl(ip)p*(y, ipfj 
+ fi(p)lz4p)+ S i(y)]p4y,p) 2 y 

(25) 

where v^ 1 — 2 d+1 n d / 2 T(d/2), derivatives of functions 
of a single argument are denoted by primes and partial 
derivatives are denoted by superscripts in parentheses; y 
is the square of the dimensionless momentum, s(y) is the 
(dimensionless) cutoff function defined from Rk(q 2 ) = 
Z k k 2 s(q 2 /k 2 ), and p*(y,ip) = [yz*{ip) + s(y) + <(<£>)] _1 
is the (dimensionless) "propagator", i.e. the 1-copy, 



2-point Green's function. Finally, dt is an operator 
acting only on the cutoff function s(y) (appearing ex- 
plicitly or through the dimensionless propagator) with 
d t s(y) = (2 — rj)s(y) — 2ys'(y). (Choices of appropriate 
functional forms for s(y) are discussed ir>i£.) In deriving 
the above equation, we have used the fact that fj = r\ and 
5* t o(ip) = z* (</?), which are properties of the dimensional- 
reduction fixed point. 

An equation for the fixed-point function is also 

derived by inserting the expansion in powers of Sip of 
(5* [see Eq. I|24p] in the corresponding beta function in 
Eq. (f23|) . The algebra is straightforward but cumbersome 
and is not worth presenting here. From the knowledge 
of u'*{ip) and z*(ip), which are obtained from two cou- 
pled equations (see RefJ^), we first solve the equation 
for S* t 2((p) and then use the input to solve Eq. (1251) . All 
partial differential equations are numerically integrated 
on a grid by discrctizing the field ip. The correspond- 
ing "cuspless" fixed point only exists above a dimension 
don — 5.15. Note also that at the upper critical dimen- 
sion, d uc — 6, one can analytically determine the solu- 
tion of the above equations: as expected the fixed point 
is Gaussian and the eigenvalue A = 1. 

The result for A versus spatial dimension d is displayed 
in Fig. „ Within numerical accuracy, the dimension d cusp 
at which A — ► is indistinguishable from do r, namely 
d cusp — doR — 5.15 (see footnote 2). We have used 
these values of A to construct the curve for the fractal 
dimension of the largest typical avalanches in Fig. 2. 
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FIG. 3: Variation with spatial dimension d of the eigen- 
value A associated with a "cuspy" perturbation around the 
"cuspless" fixed point corresponding to dimensional reduc- 
tion. Below d CU3 p — don — 5.15, only "cuspy" fixed points 
are possible. 
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61 In principle, the size of an avalanche, if defined in terms 
of the total number of spins that flip in an Ising model or 
of an equivalent geometric quantity in a field theory, could 
be different from the associated change in magnetization. 
This is for instance what is predicted for spin glasses by 
the droplet theory^. However, for the systems under study, 
one expects and proves in several cases that the size and 
the magnetization are essentially the same and scale in the 
same way with the (linear) spatial extent of the avalanche. 

62 For the equilibrium case, this is simply done through a Leg- 
endre transformation that allows one to go from the gener- 
ating functional depending on the external sources to the 
one depending on the local magnetization (the Gibbs free- 
energy in the language of magnetic systems) or the classical 
field (the effective action in the language of quantum field 



10 



theory). This can be generalized to the out-of-equilibrium 
case but requires the introduction of auxiliary fields. 
We have discussed the case where the relevant variable or 
field involved in the avalanches has a single component. 
Cusps and avalanches are also present in the case of multi- 
component fields (see e.g., Ref.— ). However, as the analysis 
is more complex, we mostly restrict ourselves here to the 
consideration of systems with a single-component field. 
In a finite-size system it may be important to appropri- 
ately sort out the various classes of large, system-spanning 
avalanches: the relevant ones for the critical scaling have 
been called "spanning critical" in the context of the out- 
of-equilibrium RFIM 3 - 3 - and d/ is then the associated ex- 
ponent. 

Weaker nonanalyticities, i.e. subcusps, could still be 



present but do not affect dimensional reductio n 13 ' 14 ). 
It has also been suggested that the critical properties of 
the RFIM in and out of equilibrium are in the same uni- 
versality class in all dimension a 25 ' 31 ' 33 
This reasoning does not apply to the elastic manifolds in 
a random environment. In this case, d is the dimension of 
space on which the displacement field describing the loca- 
tion of the manifold is defined. The size of an avalanche, 
which represents the total amount by which the manifold 
moves, therefore involves both the "lateral" dimension (by 
an amount ~ L d for a system of lateral linear size L) and 
the change in the displacement field itself (~ Z, with the 
roughness exponent £ < 1). In this case, dj = d + ( < d+1. 



